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Riddled-like Basin in Two-Dimensional Map for Bouncing Motion of an 
Inelastic Particle on a Vibrating Board 
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Motivated by bouncing motion of an inelastic particle on a vibrating board, a simple 
two-dimensional map is constructed and its behavior is studied numerically. In addition 
to the typical route to chaos through a periodic doubling bifurcation, we found peculiar 
behavior in the parameter region where two stable periodic attractors coexist. A typical 
orbit in the region goes through chaotic motion for an extended transient period before it 
converges into one of the two periodic attractors. The basin structure in this parameter 
region is almost riddling and the fractal dimension of the basin boundary is close to 
two, i.e., the dimension of the phase space. 
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Particle motion accelerated by a vibrating boundary has been studied for some time as an 
example of simple physical systems that show chaotic behavior 1, 2 ) since Fermi proposed a novel 
mechanism of acceleration of a charged particle through an oscillating field (Fermi acceleration). 3 ) 
A variant of the system is a bouncing inelastic particle on a vibrating board, and this has also 
been studied, particularly in connection with granular systems. It has been demonstrated ex- 
perimentally 4 ~ 6 ) and theoretically 7-10 ) that even a single particle bouncing on a board shows a 
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variety of behaviors such as mode-locked periodic motion, periodic doubling, chaotic behavior and 
quasi-periodic motion. 
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In this letter, based on the observation of dynamics on a single inelastic particle motion under 

•i-H , 

gravity driven by a vibrating board, we construct a simple two-dimensional map, and demonstrate 
numerically that the map shows peculiar behavior, i.e., a very long transient chaos and a riddled-like 
basin structure in the parameter region where two periodic attractors coexist. 

The original system we look at in constructing our map is an inelastic particle under the 
gravity; the particle bounces vertically on a vibrating board with a restitution constant r. Then, 
the particle velocity just before the n'th bounce u n and the velocity just after the bounce v n are 
related by 

Vn ~ V(tn) _ m 
U n ~ V(t n ) ' 1 ' 

where t n is the time at the n'th bounce and V(t) is the velocity of the board at time t (Fig.l). 
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If the motion of the board V(t) is given, we can determine, in principle, the time t n +\ and the 
particle velocity u n +\ of the n + l'th bounce from t n and v n , and we should be able to obtain a 
map that describes the particle motion by using this model. It is, however, impossible in general to 
express the map in an explicit form 10 ) because the board is moving. Here, we employ the following 
simplification that allows us to construct the explicit map by which we can calculate (v n+ i,t n+ i) 
from (v n ,t n ). 

First, we ignore the difference in the position of the board at each bounce; then we have 

2 

U n +1 = ~V n , Wl = tn + ~Vn, (2) 

where g denotes the gravitational acceleration. This simplification should be justified when the am- 
plitude of the vibration is small compared with the bouncing height (high bounce approximation). 
Under this simplification, the velocity just after the bounce must be positive such that t n +i > t n . 
To ensure this, we must assume that the vibrating board always moves upward, i.e., V(t) > 0. This 
appears to contradict the first simplification, but we may consider that the board shifts downward 
discontinuously at each period. In the present work, we employ 

V(t) = A(smujt + 1), (3) 

then we have 



v n +i = rv n + (r + l)A(smujtn+i + 1) 
2 

tn+l = t n + ~gV n . 

After rescaling the variables by the transformation, 



(4) 
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Av n , t n -> — t n , U -> — U, (5) 
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we have the two-dimensional map for v n and t n 




rv n + (1 + r)(sincjt n+ i + 1) 
t n + 2v 

n 



(6) 



with the two parameters r and uj. This was found to be equivalent to the one studied by Everson, 7 ) 
but with different parameterization. 

This map shows a typical route to chaos via periodic doubling from the mode-locking bounce 
with the same period with the external vibration as we change the parameter u> from 0.9 to 1.8 
by fixing r = 0.39 (Fig.2a). This periodic doubling and chaotic behavior have previously been 
reported for the inelastic bounce system in numerical 7 ' 8 ) and also experimental 4-6 ) studies. The 
strange attractor for the chaotic motion at w = 1.54 and r = 0.39 is shown in Fig.2b, which 
indicates that the attractor has a fractal structure with a Housdorff dimension larger than 1 in the 
two-dimensional phase space. 
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Figure 2(b) also shows a period one attractor and its basin; this period one motion is locked 
to half the vibration frequency. Upon changing the parameters further, this basin moves towards 
the chaos attractor and causes sudden destruction of the chaos through a boundary crisis when 
the basin touches the chaos attractor, then the chaos disappears as a true attractor, but its vestige 
remains as a transient one. 

After the destruction of the chaos, the period one attractor bifurcates into period two, while 
a new period ten attractor is stabilized out of the transient chaos and coexists with this period 
two attractor. In Fig.3, both attractors are shown in the phase space as well as the "temporal 
attractor" of the transient chaos which remains after the crisis. As is seen in this figure, the period 
ten attractor is embedded in the transient attractor of the chaos. In this parameter region, where 
the period two and the period ten orbits coexist, the system shows peculiar behavior, as we describe 
below. 

First, the system goes through an extremely long transient period before the trajectory con- 
verges on one of the two periodic orbits, if the initial state is chosen outside the direct basin of 
the period two attractor. The length of the transient easily exceeds one thousand iterations. The 
reason for this unusually long transient is that the region of direct attraction around the period 
ten attractor is extremely small. 

Second, the basin structures of these attractor are scattered and mixed each other almost 
completely (Fig.4). To characterize the structure of the basins, the uncertainty exponent, or final 
state sensitivity, a is introduced as the exponent for the probability /(e) that two points separated 
from each other by the distance e in the phase space are in different basins: 

m ~ t a . (7) 

The small value of a implies that the two basins are mixed each other rather wildly. The probability 
/(e) is plotted in Fig.5 for u = 1.5852 and r = 0.39, and the estimated value of a is 6.2 x 1CT 4 . 
This means the fractal dimension df of the basin boundary is very close to 2 since the uncertainty 
exponent a is related to df by 

df = d — a, (8) 

where d is the spatial dimension of the phase space. 2 ) 

The typical appearance of the basin structure of the period ten attractor in Fig.4 strongly 
resembles the riddled basin, i.e., a basin that is riddled with holes. 11 ' 12 ) Mathematically, the riddled 
basin is characterized by the following two features: (i) it has nonzero measure, (ii) any open 
neighbor of a point in the basin contains a point that does not belong to the basin. Such a basin 
has been found for a chaos attractor that is confined within a subspace of the whole phase space 
due to the symmetry of the system. According to this definition, a riddled basin is not an open 
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set, therefore the present case cannot be a riddled basin; the basins for both the period two and 
the period ten attractors are open because the attractors have a finite neighbor around them that 
is contained within the basin, therefore whole basins can be constructed by the union of infinite 
pre-images of these neighbors. 

On the other hand, the basin shown in Fig.4 cannot be distinguished from a riddled basin by 
almost any direct measure that can be calculated from the geometry of the basin. Any numerical 
evidence that the present basin is an open set is wiped away by the fragmentation of the pre-images 
during the extremely long lifetime of the transient chaos. This situation, however, does not seem 
to be unusual for the periodic window in the chaos; in such cases, the unstable periodic orbit 
embedded in the chaos attractor is stabilized extremely weakly, as in the present case. 

In summary, we have examined a simple two-dimensional map that describes an inelastic ball 
bouncing on a vibrated board, and found that the system shows a riddled-like basin structure. 
The basin should not be a riddled basin based on the mathematical definition, but its uncertainty 
exponent is very close to 2 and it almost cannot be distinguished from a riddled basin. This is due 
to the very long lifetime of the transient chaos. This situation, however, seems to be generic when 
the unstable periodic orbit in the chaos attractor is stabilized. Note that the mechanism does not 
require any symmetry in the system in contrast with the case of the riddled basin. 

The authors are grateful to Professor Kaneko for his useful comments. 
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Fig. 1. An inelastic particle bouncing on a vibrating board. 
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Fig. 2. (a) Bifurcation diagram upon changing lj with r — 0.39. (b) The chaos attractor (dots) and the coexisting 
period one attractor (open circle) with its basin (filled area) for u> — 1.54 and r = 0.39. 
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Fig. 3. The period ten attractor (solid circles) and the period one attractor (open circles) with the temporal attractor 
of the transient chaos (dots) for w — 1.5852 and r — 0.39 The orbit of the period ten attracter consists of five sets 
of two points, but the two points within each set cannot be distinguished in the figure because these two points 
are very close to each other. 
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Fig. 4. The basin of the period ten attractor (dots) for u> = 1.5852 and r — 0.39. The attractors of the period ten 
(solid circles) and the period two (open circles) orbits are also indicated. 
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Fig. 5. The logarithmic plot of the probability /(e) that two points separated by e belong to different basins, for 
u> = 1.5852 and r — 0.39. The dashed line is the line with the uncertainty exponent a = 6.2 x 1CT 4 . 



